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We present an analysis of the power fluctuations in the statistical steady state of a passively 
mode locked laser. We use statistical light-mode theory to map this problem to that of fluctua- 
tions in a reference equilibrium statistical physics problem, and in this way study the fluctuations 
non-perturbatively. The power fluctuations, being non-critical, are Gaussian and proportional in 
amplitude to the inverse square root of the number of degrees of freedom. We calculate explicit 
analytic expressions for the covariance matrix of the overall, pulse and cw power variables, pro- 
viding complete information on the single-time power distribution in the laser, and derive a set of 
fluctuation-dissipation relations between them and the susceptibilities of the steady-state quantities. 



I. INTRODUCTION 

Passive mode locking is a ubiquitous tool for creating 
ultra-short pulses, whose duration can be as short as a 
few femtoseconds. It is achieved by the simple device of 
placing a saturable absorber in a multi-mode laser cavity, 
that effectively amplifies optical waveform regions of high 
instantaneous power, thereby rendering cw waveforms 
modulationally unstable [IlUj. However, in contrast with 
dynamical model of passive mode locking, in experimen- 
tal systems with a fast saturable absorber there exists 
a threshold for self-starting of mode locking 3,6. The 
decohering noise sources present in the laser cavity are 
an inevitable impediment to the dynamical process of 
pulse buildup |4j O [7j |8]. Therefore, the self-starting 
problem could only be resolved when a statistical theory 
of interacting laser light modes subject to noise was put 
forward in 9 . In statistical light-mode dynamics (SLD) 
the optical waveform is treated as a random function, 
and the mode locking phenomenology is recovered as the 
"thermodynamics" of the SLD system. It can then be 
shown [Till E] that the noise injects entropy into the 
cavity, and that the entropy of the cw state is higher 
than that of the mode locked state. The onset of passive 
mode locking therefore obtains the significance of a first 
order phase transition between the disordered cw phase 
and the ordered pulse phase. SLD has been further ap- 
plied to predict and study several additional effects in 
passively mode locked lasers, including critical phenom- 
ena in light [12], multi-pulse phase transitions [I3J Hi] , 
noise-activated mode locking [15j . and mode locking in 
random lasers |16j . 

At the same time, noisy models of passive mode lock- 
ing were often considered, not for the purpose of under- 
standing the mode locking transition, but to study the 
temporal fluctuations of the cavity waveform [20 . The 
presence of cavity noise perturbs the output of a laser 
from a perfectly periodic train of pulses, degrading the 
pulse-to-pulse coherence, temporal periodicity, and the 
pulse shape. In applications these fluctuations are usu- 
ally obstructive, and one of the main goals of the study 
of fluctuations has been to understand how the fluctu- 
ations can be reduced. Recently the importance of un- 



derstanding and controlling fluctuations has significantly 
increased, because the pulse fluctuations are the limiting 
factor in the accuracy of the ultra-precise clocks based 
on frequency combs [TTJ (THJ QI5] • 

The standard theoretical approach to pulse fluctua- 
tions is based on perturbation theory, an apparently 
reasonable assumption, since the noise power entering 
the cavity is much smaller than the coherent power 
[2D1 [211 [53 [531 [21]. However, this argument is valid only 
for a single pass (or few passes) through the amplifier. 
In a laser that is based on a dynamical equilibrium and 
feedback, the noise power accumulates PQ, and can no 
longer be considered small. In fact a theory of fluctu- 
ations based entirely on perturbation theory can only 
capture transients. This difficulty has been recently rec- 
ognized by [2_3 [23] who demonstrated that the analysis 
of fluctuations must also take into consideration the gain 
dynamics. Therefore, a consistent theory of steady-state 
pulse fluctuations must include cw power variables in a 
statistical framework like SLD. 

The analysis of pulse fluctuations in SLD is quite simi- 
lar to that of statistical fluctuations in equilibrium statis- 
tical mechanics [28] . When there is no broken symmetry, 
thermodynamic quantities fluctuate around the equilib- 
rium values. The distribution of the fluctuations is Gaus- 
sian, and their strength tends to zero in the thermody- 
namic limit being inversely proportional to the square 
root of the number of degrees of freedom. Thermody- 
namic quantities associated with broken symmetries un- 
dergo diffusion whose rate is also inversely proportional 
to the number of degrees of freedom. There exist propor- 
tionality relationships (fluctuation-dissipation relations) 
between correlation functions of the fluctuating quanti- 
ties and thermodynamic susceptibility and transport co- 
efficients. The thermodynamic quantities associated with 
a single pulse in a passively mode locked laser are the 
pulse power, frequency, timing and phase, of which the 
last two are symmetry breaking quantities. 

In this paper we focus on the power fluctuations, since 
it is the pulse power that is directly coupled to the noise 
power and entropy through the amplifier gain [TJ 1251 126] , 
and the power fluctuations exhibit therefore the rich- 
est behavior. The power fluctuations are studied in the 
framework of a simplified model of the mode locking dy- 
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namics, where the cavity is divided into N intervals with 
the width of a pulse, N being the number of active modes 
in the cavity, and each interval is represented by a single 
complex degree of freedom. The model has been intro- 
duced and derived from the Haus master equation in |10j . 
and it has been demonstrated to yield results that closely 
approximate those obtained from the master equation 
on one hand [llj . and have large experimentally predic- 
tive content on the other hand [T3]. The model, termed 
here the coarse-grained model of passive mode locking, 
is defined in section [TTJ its thermodynamic properties are 
reviewed, and the thermodynamic power observables are 
defined and listed. The expectation values of the thermo- 
dynamic observables define the steady state of the laser, 
characterized by the overall power, the pulse power, and 
the cw power. The values of the thermodynamic quanti- 
ties are calculable from an exact mean field theory, and 
depend on the single dimensionless parameter 7 = 
where a is the absorber nonlinearity, P is the overall 
cavity power, and T is the noise power injection rate. In 
particular the pulse power is yP, where y is a function 
of 7 only. The role of 7 in SLD is analogous to inverse 
temperature in statistical mechanics, in that the order- 
ing transition occurs, and the pulse power increases, as 
7 increases [TO]. 



Since the fluctuations are Gaussian, they are com- 
pletely defined in terms of the covariance matrix of the 
thermodynamic observables, which is calculated explic- 
itly in Sec. |III| the main section of the paper. It is 
shown that, in addition to 7, the fluctuations depend on 
a second dimensionless parameter 77 = P\g'(P)\, where 
g(P) is the saturated overall net gain. 77 can be thought 
of as the gain elasticity, as it determines the restoring 
force of the gain saturation mechanism in establishing 
the steady state of the laser. If r\ becomes very large, 
the fluctuations in the overall power tend to zero, but 
there still remain pulse power fluctuations, as power is 
randomly passed back and forth between the pulse and 
the cw degrees of freedom. The covariance matrix is first 
derived directly from the invariant measure. There are 
overall seven independent elements in the covariance ma- 
trix, and the results are summarized in table|l](for entries 
depending at most on a single cw power observable) , and 
Eqs. (16-191, and convey full information on the single 



time statistical properties of the fluctuations. Next, a set 
of fluctuation-dissipation relations is derived, which are 
used to express several of the covariance matrix elements 
as derivatives of steady state thermodynamic quantities, 
thereby providing an independent verification of their 
values. Finally, the covariance matrix is also calculated 
in the much simpler case, where the laser is in the disor- 
dered phase, i.e. when it is operating in cw; there it is 
shown that the fluctuations do not depend on 7, so that 
they are unaffected by the saturable absorber. The last 
section of the paper presents our conclusions and outlook. 



II. THE COARSE-GRAINED MODEL OF 
PASSIVE MODE LOCKING 

The effective number of degrees of freedom in a mul- 
timode laser cavity is determined by the bandwidth of 
the gain medium, which determines the pulse width r p . 
Thus, a cavity of round-trip time tr has N — tr/t p inde- 
pendent degrees of freedom, which can be thought of as 
either N Fourier modes, or, in real space, as N indepen- 
dent complex field amplitudes ?/>„, each one representing 
the average value of the optical field an in interval of 
length T p . 

In this picture statistical light mode dynamics is de- 
scribed by N coupled ordinary differential equations 

mm 



d t ip n = a\4> n \ 2 ip + gip n + r/ n 



(1) 



where a is the coefficient of saturable absorption, g is the 
overall net gain, and the rj^s are independent (complex) 
centered Gaussian white noise processes with covariance 
functions (T}*{t)rj m (t')) = 2TS(t - t'). 

The equations ([!]) have to be supplemented by an equa- 
tion that expresses the gain coefficient in terms of the 
field variables using the physical gain saturation mecha- 
nism. In the case of slow gain saturation considered here, 
the gain coefficient can be assumed to be a function of the 
overall power P = jr IVvJ 2 - In many studies the gain 
saturation function is modelled by g(P) — 1+ p"/p (for a 
constant P s ) [HE]- However, it is known [2j [2T] that the 
mode locking dynamics with a |t/'| 2 '0 nonlinearity is glob- 
ally unstable with this kind of gain saturation function. 
Since, furthermore, the steady state properties of the op- 
tical waveform depend only on the local behavior of g{P) 
near the operational power P (defined precisely below), 
we will keep an arbitrary gain saturation function g(P), 
that is only assumed to have the necessary properties for 
a stable mode locking state to exist; in this way we gain 
in generality without complicating the analysis. 

The equations of motion (JT|) can be written in a vari- 
ational form [TO] 



dH 

dip* 



where H is the Lyapunov functional defined by 



2N 



£>„| 4 + NTu(P) 



(2) 



(3) 



and u(P) is the dimensionless gain potential, defined by 
-NTu'(P) = g(P). As shown in [TO], the N scaling of 
the gain potential is necessary to obtain a good thermo- 
dynamic limit as N — » 00. 

In this paper we study the steady state statistical prop- 
erties of the SLD coarse-grained model. As a white-noise 
forced variational system, the system ([!} reaches a sta- 
tistical steady state with a Gibbs-like distribution 



PbP] 



-HW/T 



(4) 
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where Z is the partition function 

Z = I [<kl>]e- H W' T 



(5) 



Steady-state expectation values of field functionals A[ip] 
are given as usual by 



(A) = / [diPWW) • 



(6) 



This is the starting point for the calculation of the mean 
values of macroscopic observables and their (co)variance. 

In previous studies [pj EH H2J H3] the thermodynamics 
of this system, i.e. the behavior of macroscopic observ- 
ables when N is large, has been studied in detail, and the 
main results are as follows: The overall power P is self- 
averaging and reaches a value P that is determined by a 
balance of the saturable absorption nonlinearity, noise 
entropy, and gain saturation. Once P is determined, 
u(P) plays no further part in the determination of the 
thermodynamics. Moreover, P itself appears in the ther- 
modynamics, except as an overall power scale, only in 
the dimensionless combination 7 = that plays here 
a role similar to inverse temperature in equilibrium sta- 
tistical physics. The thermodynamic phase diagram con- 
sists of two phases, a mode locked phase with a single 
pulse whose power Py, | < y < 1, is also self-averaging 
in the thermodynamic limit, and a cw phase where the 
intracavity power is distributed evenly between all the 
degrees of freedom. 

The thermodynamics is solvable using an exact mean 
field theory. The mean-field free energy (Landau func- 
tion) is 

f(P, y) = ^y 2 + log(P(l - y)) - u(P) . (7) 

The value of y for a given overall power is determined by 
the equation d y f(P,y) = whose solution is 



y 
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(7 > 4), 



(8) 



that does not depend directly on u and depends on P only 
through the combination 7 as explained above. There- 
fore, if only the overall power P is fixed and the pulse 
power is let to assume its equilibrium value, the free en- 
ergy becomes 



#P) = f(P,y(P)) 



(9) 



The steady state overall power is now determined by the 
equation dp(f) = that with Eq. ^ gives 



Pu'(P) 



(10) 



This equation also determines also the overall net gain g, 
since v! is proportional to it. The gain balance require- 
ment is an alternative method for reaching (10 1, based 
on dynamics |25j . 



The values of P and y derived above describe a mode 
locked state that exists for every 7 > 4, whose free energy 
is F = f(P,y). In addition, for any 7 there exists a cw 
state that can be viewed formally as a state with zero- 
power pulse, so that the mean- field free energy in the cw 
phase is 4> CW (P) = f(P, 0), and, the overall power in the 
cw state is determined by the condition dp<f> cw (P cw ) = 0, 
that is, 



P u'(P ) = 1 

± cw ^ V cw/ ^ 



(11) 



Of the two states, the statistically stable is the one with 
the lower free energy, and the other is metastable, al- 
though in general its lifetime is very large [15] . The two 
states exchange stability when 7 = 7*, above which the 
mode locked phase becomes stable. 7* is always greater 
than 4, but its precise value depends on the form of the 
gain saturation function u. For very strong gain satura- 
tion, i.e. when rj tends to infinity, 7* approaches a value 
close to 4.91 [10]. 

The overall, pulse and cw power variables are strictly 
equal to the thermodynamic quantities shown above only 
in the limit N — 00. When N is large but finite, their 
actual values deviate slightly from the thermodynamics, 
and fluctuate in time. These fluctuations are observed 
as a small randomness in the macroscopic observables. 
Since the mode system is not critical, the fluctuation 
statistics are well-described by a (multi-variate) Gaussian 
distribution [2H] . The nature of the coarse-grained model 
and the equilibrium-like approach are well suited to study 
the single-time power fluctuations. Since a macroscopic 
share of the overall power resides in the cw background 
even in the mode locked phase, the macroscopic observ- 
ables include, in addition to P and the pulse power Y , 
partial continuum power variables X a — |^n| 2 ) 

interval 

where the sum is taken on an interval containing aN de- 
grees of freedom, < a < 1 . In the interest of examining 
continuum cross-correlations we will also define comple- 
mentary continuum power variables X' b defined similarly 
to the X a variables on mutually exclusive intervals of 
length bN. 

The expectation values of all these macroscopic vari- 
ables reach finite positive limit 

(P) = P, (Y) = Py , (X a ) = aP(l - y) , (12) 

as N — > 00, while their (co)variance are 0(1/N), so that 
the relative size of the fluctuations is iV" 1 / 2 as expected. 
Our main result is the derivation of covariance matrix of 
these observables, to which we now turn. 



III. FLUCTUATION STATISTICS OF THE 
POWER OBSERVABLES 

A. The covariance matrix in the mode locked phase 

In this subsection it is assumed that the system is in 
the mode locked phase, that is, there exists a single de- 
gree of freedom whose power is close to NPy, where P 
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and y are the overall power and relative pulse power in 

the thermodynamic limit, as defined above. This aS- 
pS 

sumption implies that 7 = g ^ r - > 4. When 7 < 7* 
this state is metastable, but except for parameters very 
close to the boundary of the pulse stability region, the 
metastable state is very long-lived, and the fluctuation 
statistics are well-defined. 

The elements of the covariance matrix are the sec- 
ond order cumulants of the macroscopic observables, i.e., 
combinations of expectation values (with respect to the 
invariant measure p[ip]) of the type ((AB)) = (AB) — 
(A) (B) where A and B are two, not necessarily distinct, 
of the thermodynamic observables P, Y, X a and X' b de- 
fined above. Using relations between these observables, 
and the phase symmetry of the cw degree of freedom (see 
below), all such cumulants can be expressed in terms of 
these containing at most a single factor of X a . Expecta- 
tion values linear in X a are equal to Na times these ex- 
pectation values with X a replaced by x = |?/v!,| 2 , for n an 
arbitrary cw degree of freedom, so that we need to eval- 
uate expectation values of the form (P n Y m x k '), where 
n, m, k are nonnegative integers whose sum is at most 2, 
and k < 1. As shown in the appendix, such expectation 
values can be expressed (P n Y m x k ) = I n mk/Iooo, where 

I nmk = J dPP n+m+k J dyy m (l-y) k e N K p ^ , (13) 

where f(P, y) is the free energy defined above. The large 
N asymptotics of these integrals can be derived with the 
method of steepest decent (see appendix) and it has the 
general form 



Inmk ~ e N ^ p ^P n+m+k y m (l - y) fe ( 



1 



N 



(14) 

where P and y and / are the thermodynamic quantities 
defined above. Here ~ signifies asymptotic equality up 
to an O(l) prefactor that cancels in the calculations of 
expectation values. 

In calculating the cumulants, the leading 0(1) terms 



in I„ 



cancel, so that 



((P n Y rn x k )) = 

= F" +m+fc t/ m (l-y) fe (linear combination of J's) . 

(15) 

The technical but straightforward calculation of the func- 
tions Inmk (P, y) is outlined in the appendix, and we pro- 
ceed to state the results. 

It is possible to divide the power fluctuations into two 
parts, fluctuations that occur with a given overall power 
as a result of the random redistribution of this power 
between the various degrees of freedom, and those that 
are a consequence of the fluctuations in the overall power. 
These two contributions correspond to the 0(N~ 1 ) terms 
in the two integrations in ( 13 1 . The first contribution 



that is generated by the y integration is inversely pro- 
portional to \dyf\ = J 2 (2y — l)y, and is independent of 



the gain saturation function. It represents the dominant 
term in the fluctuations when gain saturation is so strong 
that overall power fluctuations are suppressed, that is, 
when 77 — ► 00. This term is absent in covariances that in- 
volve P. The second contribution, generated by the P in- 
tegration is inversely proportional to |9p</>| _1 = 77- 



l-2y- 



Using this information we can calculate the covariance 
matrix elements involving at most one factor of X. The 
values of the five independent matrix elements are listed 
111 table |H 

The qualitative dependence of the power fluctuations 
on the two dimensionless parameters is demonstrated in 
Fig. [T] that shows the normalized value of ((^ 2 )) as a 
function of 77 and 7. It is evident that, as expected, the 
typical fluctuations are proportional to with an 0(1) 
coefficient, when the parameters are such that the pulse 
is stable. The power fluctuations grow when the parame- 
ters approach the boundary of the of the region of stabil- 
ity, where the amplitude of fluctuations diverges. For a 
fixed value of 77, this can happen in two manners: When 
7 is too small, the pulse becomes unstable with respect 
to noise buildup, and breaks down with the laser revert- 
ing to the cw state, and when 7 is too large, the gain 
saturation is too weak to prevent the saturable-absorber 
driven growth of a pulse perturbation. Thus, for r/ > 8 
the fluctuations first decrease as function of increasing 7, 
and then start increasing as the upper stability boundary 
is approached. For fixed 7 on the other hand, the fluctu- 
ations decrease monotonically for increasing 77, reaching 
a finite value as i] — > 00. 

We now consider the remaining elements of the covari- 
ance matrix elements, containing two cw factors. These 
are 

= Jf «* 2 » + « 2 M> (16) 

((X a X' h ))=ab((xx>)) (17) 

where x = \ip n \ 2 , x' — \ip m \ , for arbitrary cw degrees of 
freedom n ^ m. Since the x variables are microscopic, 
the variance ((x 2 )) is O(l) rather than 0(N~ 1 ) and con- 
tributes significantly in Eq. (16 1. Its statistics can be 



deduced from the fact that the VVi's are (approximately) 
Gaussian random variables, and that (ip n ) — for cw 
degrees of freedom because of the phase symmetry of 



the equations of motion ([Tj), so that (x 2 ) = (iV'nl 4 ) 
2 (|^„| 2 ) 2 = 2(P(1 - y)) 2 , and therefore 



((x 2 )) = (P(l-y)f 



(18) 



The covariance ((xx 1 )} is actually linearly dependent on 
the covariances already calculated, as can be observed 
by calculating the cumulants of both sides the identity 
P - Y = Xi, giving 



((xx')) 



((P 2 ))-2((PY)) + ((Y 2 )) 

2y 



1 



1 

" N 



«* 2 » 

p2 



(2y-l)y7- 7 y 2y~l)N 



(19) 
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NP- 2 {{;■)) 
P 

Y 



Y 



(2y-i) 

(2y - l)r) - -yy 



(2y- 1)77-72/ 
f 



X„ 



(2y - 1)2/77 - iv 2 



+ 



2y - 1 V (2y - l)r) - -yy -y 2 y 



2y-l\{2y- l)r) - 



+ 



i 2 v 



TABLE I: The independent elements of the covariance matrix of the thermodynamic power observables, containing at most 
one cw power factor. The dimensionless entries in the table should be multiplied by to obtain the physical values. They 
are expressed in terms of the two dimensionless parameters 7, the inverse temperature, and 77, the gain elasticity. As r\ — > 00, 
the elements in the first row tend to zero, while those in the second row have a nonzero limit. 



B. Fluctuation-dissipation relations 

In equilibrium statistical physics, susceptibilities are 
linearly related to fluctuation covariances and transport 
coefficients are linearly related to diffusion coefficients 
[2"9"] by the property of detailed balance. These relations 
can be easily derived directly from the Gibbs distribu- 
tion. In our case, we write the Lyapunov functional as 
H[tp] = N(aR[ip] + pS[i>\), with Rty] - -3^ E„ IV>J 4 
and S[ip] = Tu(P), and let the thermodynamic quanti- 
ties and fluctuation coefficients depend parametrically on 
a and f3, with the physical quantities obtained for (3 = 1. 
In the mode locked phase the partition function is given 
by log Z = NF + 0(1), so that 



d a F 



(3. 

T 
<S> 

T 



(Py) 2 

2T 
-u(P) 



(20) 
(21) 



The left-hand equality in Eqs. (20-21 1 follows by differ- 



entiation of the partition function under the integral sign, 
and the definition ^ of the expectation values, and the 
right-hand equality from direct differentiation of Eq. ^ . 
The results obtained for (R) and (S) are actually special 
cases of the fact that (k(P,Y,X a )) = k(P , Py,aP(l-y)) 
for any function k of the thermodynamic observables in 
the thermodynamic limit, that holds since they are self- 
averaging. Note that R ~ -V in the thermodynamic 
limit, since the cw terms in the sum give an O(jj) con- 
tribution. 

Taking another derivative of the free energy then gives 
the second order cumulants, 



daFmi 
d a dpF m \ 



dgFmi — 



((R 2 )) _ d a (R) 



rp2 

({RSI 
y2 



T 

d a (S) 
T 

T 



(22) 
(23) 
(24) 



The right-hand equalities in Eqs. ( 22 24 ) are obtained by 



replacing the first derivative by its value from Eqs. (20 
21 1. In the mixed derivative case there are two possible 
replacements that give the Maxwell relation d a (S) = 
dp(R). 



Once more we can use the self-averaging of P and Y 
simplify the cumulants; in this case it allows us to ex- 
press them in terms of the elements of the covariance 
matrix. For general functions k\ and ki self-averaging 
implies that 



«fc 1 (P, Y)k 2 (P, Y))) = d P k 1 dpk 2 <(P 2 )) 

+ (dpk 1 d Y k 2 +d Y k 1 d P k 2 ) ((PY))+dyk 1 dyh 2 {(Y 2 )) 



(25) 



where the function on the right-hand-side are evaluated 
as usual at P = P, Y — Py. It follows that 



((R 2 )) = (Pyf((Y 2 )) , 

((RS)) = PyTu'(P) 2 ((P 2 )) , 
{{S 2 )) = T 2 u'(P) 2 «P 2 )) . 



(26) 
(27) 
(28) 



Thus, Eqs. (22 24 1 offer an alternative method to cal- 



culate three elements of the covariance matrix derived 
in Sec. |III A| by taking the parametric derivatives of R 
and S, d (R) = ^(Pdy + ydP), OS = Tu'(P)dP, where 
the '<9' stands for differentiation with respect to either 
parameter. The parametric derivatives of P and y are 
obtainable by short calculations starting from Eqs. ^ 
and (10), and give the susceptibilities 



dpP 
d Q P 



l- 2 y 

d P 



2y 



1 P 

P 2 y 



7 T(2y-l) 
1 /2 7 



dpP 



d a P 



p2 

T 



(29) 

(30) 
(31) 
(32) 



7 2 (2y — 1) V P 

The relevant elements of the covariance matrix are repro- 



duced when Eqs. ( 26 32 ) are substituted in Eqs. (22-241 



C. Fluctuations in the cw phase 

For completeness we also present the power fluctua- 
tions in the cw phase, where power is evenly divided 
(on average) between all the laser degrees of freedom. 
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FIG. 1: A contour plot of the pulse power fluctuations vari- 
ance ((Y 2 )) as a function of the dimensionless parameters 7 
(inverse temperature) and r\ (gain elasticity). Lighter shades 
represent stronger fluctuations and darker weaker fluctua- 
tions, except the dark blue region in the left part of the figure, 
that designates parameter values where the pulse is unstable. 

Contour lines represent the values (from right to left) 0.2, 0.3, 

pi 

0.5, and f, times 



This state is stationary for 7 < 7* and metastable for all 
7 > 7*, although for very large 7 it is short-lived. The 
latter case is the self-starting regime [5J [7J [TS] , and our 
results hold away from this regime. 

Because the saturable absorption is nonlinear, F cw — 
/(P C w, 0), the cw free energy, is independent of 7 (see Eq. 
(eq:fpy)), and therefore the cw thermodynamics is inde- 
pendent of 7 [TU]. This is also true for the fluctuations, 
that are therefore (up to small corrections) those of free 
randomly forced modes subject to the gain saturation 
constraining potential u. The calculation of the covari- 
ance matrix is very similar to that of the more compli- 
cated mode locked case, and we simply state the results. 

The power variance is 



•< P »c„ N \d 2 P cf> cw (P cw )\ 



(33) 



Since moments including Y are by definition equal to 



zero in the cw state Eqs. (16 19 1 together with Eq. (33) 
determine the remaining nontrivial cumulants (IX^jj 
and 



IV. CONCLUSIONS 

The main result of this paper is a complete quantita- 
tive characterization of the Gaussian fluctuations in the 



overall, pulse, and cw power of passively mode locked 
lasers. It is the first application of SLD to the study of 
steady-state fluctuations, and further demonstrates the 
power of this method in the understanding of the interac- 
tion of noise and nonlinearity, and its consequences. The 
analysis reveals deeper consequences of the close analogy 
between a mode locked laser and the equilibrium statis- 
tical mechanics of interacting mode systems. In particu- 
lar the analogy led us to identify and prove fluctuation- 
dissipation relations for the far-from-equilibrium laser 
system. 

The results were obtained in the framework of a sim- 
plified mode locking model. The coarse-grained model 
provides a good approximation for the behavior of more 
realistic models and experimental systems, but the quan- 
titative details are model-dependent. Nevertheless, sev- 
eral important conclusions can be drawn from the quali- 
tative properties of our results, that are likely to hold for 
a wider class of systems. First, we observe that the fluc- 
tuations are determined by a further dimensionless, the 
gain elasticity, in addition to the parameters that deter- 
mine the steady state (here a single parameter) . The gain 
elasticity measures the response of the amplifier gain to 
changes in the overall power, and therefore controls the 
strength of fluctuations in the overall power. As gain 
elasticity decreases, overall power fluctuations increase, 
as well as fluctuations in the pulse power and cw power. 
Secondly, we find that the overall power and the pulse 
power are always positively correlated, while cw power 
and pulse power are negatively correlated. This informa- 
tion can in principle be used to predict and correct for 
pulse fluctuations before the arrival of the pulse. On the 
other hand, cw power accumulated in different parts of 
the cavity display positive correlations. 

As in all noise driven systems, the fluctuations diverge 
as the parameters approach the boundary of the region 
of stability. Close to this boundary, fluctuations cease 
to be Gaussian; this case can also be studied within the 
SLD framework [15], but with more sophisticated meth- 
ods than those used here. In particular, the pulse is sus- 
ceptible to elimination by noise-activation transition to 
the disordered state. Another potential region of inva- 
lidity of the Gaussian distribution, is the neighborhood 
of a critical point [SB], absent in the model studied here, 
but present in more general laser systems [15] . In either 
case a full analysis of the light fluctuations has to take 
into account the quantum properties of the spontaneous 
emission noise [3D]. 

The most serious limitation of the static method used 
here to calculate fluctuation properties is that it can only 
capture single-time statistics. Many-time statistics, like 
timing and phase jitter can also be naturally studied 
within SLD, but, like in statistical mechanics, they re- 
quire consideration the kinetics of the electromagnetic 
field, in addition to its steady-state properties. This task 
is postponed to future study. 
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APPENDIX A: STEEPEST-DECENT 
EVALUATION OF COVARIANCE MATRIX 
ELEMENTS 



evaluation of cumulants 
The expectation val- 



III A 



Here we outline the 
((P n Y m x k )) defined in Sec 

ues are calculated starting with the definition Eq. (|6]) 
fixing the overall power P and first integrating over the 
N — 2 degrees other than the cw degree of freedom ip n , 
letting Px = \ip n \ 2 , and the pulse degree of freedom ip p , 
with NPy = iV'pl 2 , obtaining 

x J dxx k (P(l-y- x/N)) N , (Al) 



where the term (P(l — y — x/N)) N , the result of the N — 2 
integrations, is equal to the cw-state partition function 
up to an irrelevant 0(1) pref actor [TO]. The x integration 
then leads to Eq. ( 13 ) defining I nm k- The remaining inte- 



grals over y and P are then evaluated for large N using 
the standard and well-known steepest-descent method, 
otherwise known as the Laplace integral method [21] ■ It 
is based on approximating the integrand with its Taylor 
series near the maximum of the exponential at (P,y). 
We choose to perform the integration in two steps. The 
result of the y integration is I nm k — J dPI nm k(P), 



I ,( P ) - J 2 JL e Nf(P,V) pn+m+k ( -rn (1 _ -sk , JL [ vf%_ \ Q ( J_ \| ( M) 

where the • • • stand for several other O(jj) terms that cancel in the evaluation of the cumulants. 

The O(jj) term in Eq. (A2) is carried into the next integration, over the P variable, after which P is replaced by 
P, and a similar 0(jj) is generated giving 

j _ d 2 y (y m (l-y)% {P) d P (P n + m + k y(Pr(l~y(P)) k )\p f ^ 

W- W /(P,y)| + 2\dU(P)\ • [A6) 
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